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Abstract 

Lam’s phenomenological e-renormalization group (RNG) model is quite different from 
the other members of that group. It does not make use of the correspondence principle and 
the e-expansion procedure. In this report, we demonstrate that Lam’s e-RNG model [Phys. 
Fluids A, 4 , 1007 (1992)] is essentially the physical space version of the classical closure 
theory [Leslie and Quarini, J. Fluid Mech., 91 , 65 (1979)] in spectral space and consider the 
corresponding treatment of the eddy viscosity and energy backscatter. 


^This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASI- 19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, \A 23681-0001. 
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Introduction 


In this note, we demonstrate that Lam s e-RNG modeL is essentially the physical space 
version of the classical closure theory^ in spectral space and consider the corresponding 
treatment of the eddy viscosity and energy backscatter. 


Analysis 

The incompressible N-S equations are 

^ + V • Vv = --Vp + i/qV^v (1) 

Ot p 

where vq is the molecular viscosity, p is the density, and p is the pressure and can be 
determined from (1) using V • v = 0. The external driving force that sustains the turbulence 
and which acts in the very small wavenumber region is not included in (1) since it plays no 
part in the energy cascade process in the inertial range^. 

As in both closure and RNG theories, the velocity field is filtered into two components 

v = v"^ + v^, p = p'^+p^ (2) 


where the Fourier-transformed fields 


v<{kJ) = G{k)v,{k,t), (3) 

v>{k,t) = {l-G(k))v,{k,t). (4) 


The sharp cut-off filter of classical closure theory is exactly the same as the RNG technique 
of separating the subgrid from the resolvable scales at the cutoff wavenumber A 


In the classical closure theory 


G{k) = I 
of Leslie and 


0 if A: > A ; 

1 if k < A. 

Quarini (LQ)^, the filtered N-S equation is 


(4) 


+ [,,o + UE{k)]k^)v<{kJ) = M.0-y{k) J dpdqu|(p,L)<(q,L) -f f.(kj). (5) 
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where Mae-yik) is the standard nonlinear coupling coefficient^’^. For convenience we have 
added to both sides a wavenumber dependent turbulent eddy viscosity VE^k). which is at the 
moment unspecified. The term f{k^t) accounts for the Reynolds stress^’^, 

R/3^ = ( 6 ) 

the cross stress^’^, 

Cg^ = n^(p,<)n>(q,t) + u|(p,t)u<(q,t) (7) 

and the added eddy viscosity u^ik) : 

/c(k, t) = OE{k)k^v<{k, t) + Mcp-y{k) J dpdq [Cp-^ + Rp-y] . (8) 

In (6)-(7), |p + q| < A. It is important to realize that no random force has been inserted 
here. 

In the Lam approach to e -RNG^, one works in physical space rather than wavenumber 
space. The exact resolvable scale Navier-Stokes equations can be written 

O A 

[- - (I/O + i/t)V']v< = ~-Vp< - V • (v<v<) + (9) 

where is defined by 

gfast _ y _ (v<V< - vv) — i/jV^V"^ = V • (2v^V< — V^V^) — I/yV^V. (10) 

Note that Lam has introduced a k-independent turbulent eddy viscosity, i/j .which remains 
to be chosen. is generated by the filtering process. The term g-^“®* in physical space 

corresponds to the term f(k, t) in wavenumber space, in Eq. (8). 

The classical theory proceeds from this point by the use of certain “closure approximations”^ 
An equation for the resolvable spectral energy, E{k, t ), can readily be derived, 

- + 2i^ok^]E{k,t) = T{k,t) + T>{k,t), (11) 
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where T{k,t) is the resolvable scale energy transfer and T>{k,t)\s the energy transfer caused 
by the cross and Reynolds stresses^ which can be put into the form^’’^ 

T>(k,t) = -2Mk)k^EikJ) + Uik). (12) 

U{k). which represents the backscatter of energy from small to resolvable scales and is also 
the spectrum of the correlation function of f, is given by 

U{k)= I dpdqB[k,p,,,]E(p)Elq)G^k)ll-G{p)G{q)]. (13) 

Ja 

Vd{kA). the drain eddy viscosity, is given by 

i/d{k) = f dpdqA{k,p,q)E(q)[l - G(p)G(q)]. (14) 

Ja 

The integration domain is denoted by the expression A in which p and/or q > A. The 
explicit functional forms of A and B appearing in (13)-(14) are given in Leslie^ and LQ^. 

Instead of trying to compute using closure approximations, Lam^ simply tries to 

model its correlation function based on physical arguments. In his view, f is simply a guess 
of what should be for k ^ A in the resolvable scale Navier-Stokes equation. He noted 

that in the absence of f, the energy spectrum of the flow, computed from (5) driven by initial 
and/or boundary conditions, will have a Kolmogorov dissipation wavenumber substantially 
smaller than A. The primary role of f is to extend for the resolvable scale velocity field the 
inertial range with a guaranteed Kolmogorov scaling for k ^ A and beyond. 

The forcing function in classical closure theory arises from filtering at the small scales. 
In modeling the correlation function of f , Lam' assumes the form 

< /,(k,^)/,(k',u.') >= J-£-^r''+'-<(2r)">P.,(k)6(k+k')«(u. + u,.') (15) 

where w is frequency, S is the dissipation rate, d is the dimension of the physical space, Ha is 
a constant, and P,jik) = <5,j - kykjjk'^. A multiplicative factor involving A**”' is introduced 
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to maintain dimensional consistency for arbitrary c. It is of some interest to compare Eq. 
(15) with the forcing correlation function introduced by Yakhot and Orszag (YO)^ 


< MKu:)f,{k\u;') >= ^ek-^^‘^-^i2KY+^P,,{k)6{k + k')6{u, + J), (16) 

where 0 is a known constant determined by 2DoSdl(2TrY'^^ = 1.594£’ (YO®) and Sd is the 
area of a (/-dimensional unit sphere. This form' is assumed to arise from forcing at k = 0 ; 

<ff>^6{k)£6{k + k') (17) 

with the use of Gel’fand’s ^-function representation in the limit of e — 4 and A; — > 0 

6{k) = lim(4 — e)A:^“^ for k 0 (18) 

e— ^-4 ' ^ 


To recover (16), it appears that (18) needs to be applied for ^ 0, without the (4 — e) factor. 

Lam pointed out that the forcing correlation function, Eq. (15), should peak around 
A; that its magnitude should be small for small k by an appropriate choice of i/t', and 
that its behavior for k >> A is unimportant and irrelevant for the evolution of the resolved 
modes. Most importantly, the correlation function now depends on A, while in e-RNG®“”, 
the correlation function is assumed to be “scale invariant”. The dimensionless parameter e 
in the correlation function is now available as a freely adjustable parameter, and Lam used 
it to make the “predicted value” of Kolmogorov constant acceptable. He showed that either 
e = 0 or e = 0.923 yield good results. 

The stochastic backscatter f, for isotropic homogeneous turbulence in three dimensions, 
has a k"* spectrum to lowest order in wavenumber k (e.g., Ref. 5). Specifically, 

U{k) = ^k^ for ^ 0. (19) 

where 0k,p,q{t) = ?,,(/) + + q^)] and Hk,p,q{t) is an “eddy-damping rate” of the 

third-order moments associated with the wavevectors k, p, and q. 
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Thus, Lam’s postulate (which was based on intuitive physical arguments) that U{k) is 
small for small k is consistent with classical closure theory. 

The advantage of the classical theory is that the energy equation is always satisfied and 
no restriction on the magnitude of A is imposed — so long as A is in the inertial range. On 
integrating (11) with respect to A' for 0 < A: < A, we obtain: 

^ = n-E. ( 20 ) 

ot 

where K is the integral of E{k) over the resolved wavenumbers, and S is defined by: 

£= f^T>{k)dk= l^2k^Un(k)E{k)dk. (21) 

Jo Jo 

and n, the resolved energy transfer term, is given by: 

n = l^f{k)dk. 

Jo 

The net eddy viscosity, Un{k,t), is defined^’^’®“^ as 

Vn{k) = Vd{k) - Ub{k). (22) 

and i/f){k,t), the back-scatter viscosity, is given by 

Vb{k)^U{k)l{2k'^E{k)). (23) 

From (14) and (23), one can show^° that for k in the inertial range and k <C A, the ratio of 
Vb{k) to Ud{k) is equal to j|(^/A)“'^^. Spectral large-eddy simulations (LES) of Lesieur and 
Rogallo^’^^ was based on the resolvable scale Navier- Stokes equation 

[■^ + Wo + i^n{k)]k^)v<{k,t) = Mc 0 -,{k) j j dpdqu|(p,t)f<(q,/). (24) 

Lam emphasized that £, the energy dissipation rate of the turbulent flow in question, 
must be related to the parameters of the turbulent eddies by an ad hoc postulate under his 
formulation. Lam’s choice* is 
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( 25 ) 


Sl = lim 2i/t{A) [ k^E{k)dk. 

A-^oo Jo 

The large A limiting process in (25) is needed to ensure that the dissipation rate can be 
adequately evaluated using information available from the resolved modes alone. In Lam’s 
approach, the value of A must be sufficiently large such that the dissipation function El as 
given by (25) is independent of A. In physical variables, El is defined by: 

= ■ (26) 

The Smagorinsky result for vj is recovered if El is eliminated between (26) and J^r(A) = 
CyE^L^ . In LES, the Lam requirement that A must be large enough is equivalent to 
requiring that (26), computed using data only from resolved modes, be “grid size” inde- 
pendent. In Lam’s view, an LES calculation must exhibit a Kolmogorov spectrum using 
the resolved modes such that the limiting process in (25) is respected. If it does not, then 
the calculation would have no theoretical standing. Physically, if A is sufficiently large (so 
that El is independent of A), the contribution of back scattering to the dissipation would be 
negligible. The random force f , the surrogate of the does not appear explicitly in the 

final LES model of Lam and one needs only to provide a profile of < ff > so as to introduce 
the adjustable parameter e used in computing i/x- 

Conclusion 

Thus, we find that Lam’s formulation of e-RNG^ is essentially the physical space version 
of the spectral classical closure theory^ with Vn{k) being replaced by a phenomenological 
^-independent i/t, but which now depends on arbitrary parameter e. 
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